Bound states in weakly deformed 
strips and layers 
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Abstract. We consider Dirichlet Laplacians on straight strips 
in M 2 or layers in R 3 with a weak local deformation. First 
we generalize a result of Bulla et al. to the three-dimensional 
situation showing that weakly coupled bound states exist if 
the volume change induced by the deformation is positive; we 
also derive the leading order of the weak-coupling asymptotics. 
With the knowledge of the eigenvalue analytic properties, we 
demonstrate then an alternative method which makes it possi- 
ble to evaluate the next term in the asymptotic expansion for 
both the strips and layers. It gives, in particular, a criterion for 
the bound-state existence in the critical case when the added 
volume is zero. 

1 Introduction 

Spectra of Dirichlet Laplacians in infinitely stretched regions such as a 
planar strip or a layer of a fixed width have attracted a lot of attention 
recently. Of course, the problem is trivial as long as the strip or layer is 
straight because then one can employ separation of variables. However, 
already a local perturbation such as bending, deformation, or a change 
of boundary conditions can produce a non-empty discrete spectrum. 

This effect was studied intensively in the last decade, first because 
it had applications in condensed matter physics, and also because it 
was itself an interesting mathematical problem. A particular aspect 



we will be concerned with here is the behaviour in the weak-coupling 
regime, i.e. the situation when the perturbation is gentle. 

Recall that the answer to this question depends on the type of 
the perturbation. For bend strips, e.g., one can perform the Birman- 
Schwinger analysis which yields the first term in the asymptotic ex- 
pansion for the gap between the eigenvalue and the threshold of the 
essential spectrum flDEfl . It is proportional to the fourth power of the 
bending angle and always positive, since any nontrivial (local) bending 
induces a non-empty discrete spectrum. A local switch of the bound- 
ary condition from Dirichlet to Neumann has a similar effect. Here 
the weak-coupling behaviour was determine variationally to be gov- 
erned by the fourth power of the "window width" [EVTJ and the exact 
asymptotics was derived formally in |Po|1 by a direct application of the 
technique developed in |TJ, |Ga| |. Notice that this asymptotics differs 
substantially from that corresponding to a local change in the mixed 
boundary conditions, where the Birman-Schwinger technique is ap- 
plicable and the leading term is a multiple of the square of the said 
parameter |[EK|| . Recall also that analogous results can be derived for 
layers with locally perturbed boundary conditions where, however, the 
asymptotics is exponential rather that powerlike [|EV2| . 

The present paper deals with the case of a local deformation of 
the strip or layer, which is more subtle than the bending or boundary- 
condition modification. The main difference that the effective interac- 
tion induced by a deformation can be of different signs, both attractive 
and repulsive. It is easy to see by bracketing that a bulge on a strip or 
layer does create bound states while a squeeze does not. The answer is 
less clear for more complicated deformations where the width change 
does not have a definite sign. 

The first rigorous treatment of this problem was presented in the 
work of Bulla et al ||BGRS|| dealing with a local one-sided deforma- 
tion (characterized by a function Xv ) of a straight strip of a constant 
width d. The authors found that the added volume was decisive: a 
bound state exists for small positive A if the area change Xd(v) is pos- 
itive, and in that case the ground-state eigenvalue has the following 
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weak-coupling expansion, 



E(X) 



A*K*<v) a + 0(A d ) 



where Ki = ^ is the square root of the first transverse eigenvalue.^] On 
the other hand, the discrete spectrum is empty if (v } < 0. A problem 
arises in the critical case, (v) = 0, when the areas of the outward 
and inward deformation coincide. The authors of BGRS | suggested 
that the analogy with one-dimensional Schrodinger operators by which 
bound states should exist again may be misleading due to the presence 
of the higher transverse modes. 

This suspicion was confirmed in ||EV3|| where it was shown that 
this is true only if the deformation was "smeared" enough. More specif- 
ically, the discrete spectrum is empty if 



d> —=b 



;i.2) 



provided suppt> C [—b,b]. On the other hand, a weakly bound state 
exists if 



,/||2 



< 



9 + V90 + 12tt 2 
and in that case there are positive c±, c 2 such that 

-ciA 4 < E(\) — k\ < -c 2 A 4 . 



1.3) 



;i-4) 



These results have been obtained by a variational method and they are 
certainly not optimal, because there are deformed strips which fulfill 
neither of the conditions (|1.2|), (|1.3|). 

A way to improve the above conclusions would be to compute the 
Birman-Schwinger expansion employed in [|BGRl5|] to the second order 



which becomes the leading one when the term linear in A 2 in ( |1 . 1|) 
is absent, and the asymptotics is governed by A 4 in correspondence 
with (|1.4|). This is not easy, however. The standard technique in these 



*In fact, they assumed d = 1, but it is easy to restore the strip width in their 



expression obtaining eq. (1.1) 



3 



situations is to map the strip in question onto a straight one by means 
of suitable curvilinear coordinates. In distinction to the bent-strip 
case [DE| these coordinates typically are not locally orthogonal. Hence 
the transformed Laplacian contains numerous terms which make the 
computation extremely cumbersome. 

After this introduction, let us describe the aim and the scope of 
the present paper. The aim is twofold. First we are going to consider 
an extension of the result of ||BGRS|| to the case of a locally deformed 
layer. The result is summarized in Theorem |2]4]. In particular, we 



derive a weak-coupling expansion of the ground-state eigenvalue, 



E(X) = k\ — exp 



2 (-\^( V ) + 0(X 2 ) 



7T 



;i.5) 



and show the analytical properties of the round-bracket expression 
w.r.t. A. This is done in Sec. 2; the results again say nothing about 
the behaviour in the critical case. 

Instead of attempting to proceed further by the Birman-Schwinger 
method, we demonstrate in Sec. 3 a different approach to the weak- 
coupling problem. It is based on constructing the asymptotics of a 
particular boundary value problem, and requires as a prerequisite the 
analyticity of the function E(-) itself in dimension two, and of its 
above mentioned constituent in dimension three. In the present case, 
however, these properties are guaranteed by |[BGRS|| and the results of 
Sec. 2. The methods allows us to recover the expansions ( |1 . 1| ) and ( |1.5 ) 
in a different way. What is more, we are also able to compute higher 
terms, in principle of any order. We perform the explicit computa- 
tion for the second-order terms which play role in the critical case. In 
particular, we made in this way more precise the result expressed by 
( |1.2j ) and ( |1.3| ) about the critical bound-state existence for smeared 
perturbations, and derive its analog in the deformed- layer case. 



2 Locally deformed layers 
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2.1 The curvilinear coordinates 



Let x = (x 1 , x 2 ) E 



a function v G 



(2.2) 



2 and (x, u) E Q := R 2 x (0, d) with d > 0. Given 
R 2 ) we define the mapping (A > 0) 

: -> M 3 : {(a:,u) i-> 0(a;,u) := (r\a; 2 , (1 + Xv(x))u)} (2.1) 

for A > 0, which defines our deformed layer f2 A := <j)(£l ). 

To make use of the curvilinear coordinates defined by the map- 
ping we need the metric tensor := 4>,i-4>,j of the deformed layer. 
It can be seen easily to be of the form 

(1 + X 2 v\u 2 \ 2 v t iv !2 u 2 
\ 2 V tl V,2U 2 1 + \ 2 V%U 2 

A-u i(l + Xv)u Af j2 (l + Xv)u 

where v ttl means the derivative w.r.t. and its determinant is G := 
det(G ij ) = (1 + Xv) 2 . 

In view of the inverse function theorem, the mapping <fi defining 
the layer will be diffeomorphism provided AH^-Hoo < 1, where we 
put conventionally V- := max{0, — v}. For a sign-changing v, this is 
a nontrivial restriction which is satisfied, however, when A is small 
enough. That is just the case we are interested in. 

We will also need the contravariant metric tensor, in other words 
the inverse matrix 

/ 1 \ 



Aui(l + Xv)u 
Xv 2(1 + Xv)u 
(1 + Xv) 2 



(G ij ) 



V 







At; in 



Xv 211 



1+Xv 1+Xv 

and the following contraction identities 



1+Xv 



Xv^u 
~ l+Xv 

1+A 2 |Vd| 2 m 2 
(1+At,) 2 



(2.3) 



Xv , 



XAvu 3X 2 \X7v\ 



+ 



'U 



(2.4) 



1 + Au' " J 1 + Au ' (1 + Av) 2 ' 
where conventionally summation is performed over repeated indices, 
and we denote | Vt>| 2 := v\ + v\ and Av := v t u + f,22- Another con- 
vention concerns the range of the indices, which is 1,2 for Greek 
and 1,2,3 for Latin indices. The indices are at that associated with 
the above coordinates by (1,2,3) <-> (x 1 , x 2 , u). 
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2.2 The straightening transformation 

As mentioned in the introduction the main object of our study is 
the Dirichlet Laplacian — A^, A on L 2 (Q\). If we think of a quantum 
particle living in the region Q\ with hard walls and exposed to no 
other interaction, — A^, A will be its Hamiltonian up to a multiplicative 
constant; we can get rid of the latter by setting the Planck's constant 
h = 1 and the effective mass m* = ~. Mathematically speaking, — A^, A 
is defined for an open set C M 3 as the Friedrichs extension of 
the free Laplacian with the domain C^°(fi) - cf. [ |RS| , Sec. XIII. 15]. 
Moreover, since the smooth boundary of Q\ has the segment property, 
— A^, A acts simply as ip i— > —if>jj with the Dirichlet b.c. at dQ\. 

A natural way to investigate the Hamiltonian is to introduce the 
unitary transformation U : L 2 (Q\) — > L 2 ({Iq) : {ip i— ► Uip := Giipocfy] 
and to investigate the unitarily equivalent operator 

H x := Ui-A^U- 1 = -G-^diG^G^djG- 1 * (2.5) 

with the form domain Q(H\) = Wq 2 (Q ) instead of — A^, A . As usual 
in such situations, the "straightened" region is geometrically simpler 
and the price we pay is a more complicated form of the operator ( |2.5| ) . 

To make it more explicit, put F := InG^. Commuting G~* with 
the gradient components, we cast the operator ( 2.5 ) into a form which 
has a simpler kinetic part, 

H x = -diG ij d 3 + V = -G^didj - (;'■[ <), + V , 

but contains an effective potential, 

v : (('■■' I-'.,).. + /•'/'"•'/•:, = G'-'i-'.ij + (■ {j-'. i + c" j /•;,/• :,. 



If we now employ the particular form Q2.2|) of the metric tensor to- 
gether with ( |2.3| ), (|2T4j), we can write 

2 2 l + X 2 \Vv\ 2 u 2 2 2\v±u 2Xv, 2 u 



(1 + Xv) 2 1 + Xv 1 + Xv 



Xv A Xv, 2 a , fXAvu 3A 2 |Vv|V. ... 

1 + Xv 1 + Xv V 1 + ^ (1 + 
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with 



V 



A Aw X 2 vAv 3A 2 |Vw| 



2 2(1 + Aw) 4(1 + Aw) 2 ' 

For our purpose it useful to rewrite this expression further in a form 
sorted w.r.t. to the powers of A: 



H x = -Ag° + A 



2w d 2 + 2w iiidids + 2v 2 ud 2 d 3 + v i<9i + vA 



, . . _ Aw 
+ (Aw) w<9 3 + — 

3w 2 + |Vw| 2 w 2 + 2Aw 3 o9 2vv 1 u _ _ 2vv 2 u _ _ 



(1 + Aw) 5 



1 + Aw 



1 + Aw 



w,i o vv j2 a , ('v(Av)u 3|Vw| 2 w 

+ ^^ + —^+(^- + ^-^3 

wAw 3|Vw| 2 
+ ^ r-r + 



1 + Aw (1 + Av) 5 



2(1 + Aw) 4(1 + Aw) 2 



In analogy with [ |BGRS|| , we thus get the following formula for the 
"straightened" operator, 



H X = H + \Y, KBn + A 2 A nBn 



(2.6) 



n=l 



n=4 



where each of the A^s and B n : s is a first-order differential operator 
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with compactly supported coefficients and 

A\ := 2vd 3 B x := cod 3 

A\ := Av B 2 := uo ( ud 3 + 1 



2 

A* 3 := (2ud 3 + B 3 := v tl di + v t2 d 2 

3v 2 + |Vv|V + 2Av 
(1 + \v 3 ) 



A l ■= k \ * d * B A := cud 3 



vAv r> ( o 1 

^ := -IT^ B,:=.[ud 3 + - 

A l ■= - ,^ V yL 9 B 6 := uo (ud 3 + ~ 



:i + Av) 2 ' V 4 

2w<9 3 + l 

A := — i — ; — \ — 11 5 7 := V, 1#1 + ^,2<9 2 

1 + Xv 

with u; G C^°(IR 2 ) such that u = 1 on suppf. We define a pair of 
operators C x , D : L 2 (tt ) -> £ 2 (fi ) ® C 7 by 



A„y? n = l,2,3 
AA n <^ n = 4, . . . , 7 



<p ^ (D<p) n := B n <p n = l,..., 7 
then (O) finally becomes H x = H + \C* X D. 



2.3 Weak coupling analysis 

First we note that since the our layer is deformed only locally, we have 

<7 ess (-Ag*) = (7 ess (-A£°) = [k 2 , oo) . 

This is easy to see, for instance, by using a bracketing to show that 
inf a ess (— A^, A ) = k\ - cf. [pEKf - while the opposite inclusion is 
obtained by constructing an appropriate Weyl sequence. We use the 
notation k, 2 := (^j) 2 for the eigenvalues of the transverse opera- 
tor (— d 2 ) D ; the corresponding eigenfunctions are denoted by Xji an d 
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their explicit form is 



f \ 2 ■ 
XjW = Y -i sin K n u . 

Next we define K% := XD(H Q — a 2 )~ x C^. We are interested in (posi- 
tive) eigenvalues E(X) =: a 2 of if^ below the lowest transverse mode, 
hence we choose a G [0, Ki). Our basic tool is the following classical 
result - cf. [[BGRS| , Lemma 2.1]: 

Proposition 2.1 (Birman-Schwinger principle) 



« 2 e o"disc(^A) <^ -1 e o"disc(^. 



Proof: If = -if), then <p := -X(H - a 2 )' 1 ^ is easily 

checked to satisfy H\ip = a 2 ip. Conversely, if H\(p = a 2 ip, we have 
<p G Q(H X ) C D(D), so ip := is in £ 2 (ft ) and = I 

To make use of the above equivalence, we have to analyze the 
structure of K%. Let Ro(a) := (Hq — a 2 )" 1 be the free resolvent cor- 
responding to Hq. Using the transverse- mode decomposition and the 
fact that Hq = — A R g) I\ + I2 ® (— <9f ) D , we can express the integral 
kernel of R , 

00 

Ro(x, u, x', u'; a) = Xj{ u ) r j( x > x '~i a ) Xj{ u ') 
3=1 

where rj(x, x'\ a) is the kernel of (— A r2 + k 2 — a 2 )" 1 in L 2 (R 2 ). We 
define kj(a) 2 := k 2 — a 2 . The free kernel can be expressed in terms 
of Hankel's functions - cf. [ |AGH| , Chap. 1.5] - which are related to 
Macdonald's functions by |[AS| , 9.6.4], so finally we arrive at the formula 

R (x,u,x\u';a) = — J^Xi(w) K (kj(a)\x - x'\) Xji u ') ■ 

i=i 

Now we want to split the singular part of Rq ; we write = L\ + M\ 
where L\ := \DL a C{ contains the singularity: 

L a (x,u,x',u) := - — xi{u) lnki(a)xi{u') 

Z7T 
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diverges logarithmically as a — > K\—. The regular part M A = XDM a C^ 
consists of two terms, M a = N a + R$ (a), where the operator R$ is 
defined as the projection of the resolvent on higher transverse modes 

^ oo 

Ro(x, u, x\ u'; a) := — ^ Xj(u) K^k^a^x - x'\) Xj(u'), 

71 3=2 

and the remaining term is therefore 

N a (x,u,x',u') := ^-Xi(«) (^oOiO)^ - x'\) + In k x (a) )xi(V)- 



Put w^ 1 := hxki(a). The next step in the BS method is to show the 
boundedness and the analyticity (w.r.t. w) of the regular part of K®. 
A more difficult part of this task concerns the operator containing N a 
where we have to take a different route than that used in |BG11S|| . 



First we note that while the Hilbert-Schmidt norm is suitable for 
estimating the operator N a , it fails when the latter is sandwiched be- 
tween XD and C^. More specifically, using the regularity and compact 
support of the functions involved one could transform \DN a C^ into 
an integral operator via integration by parts, but the obtained kernel 
has a singularity which is not square integrable. Hence we use instead 
the "continuous" version of the Schur-Holmgren bound. Since it seems 
to be less known than its discrete analogue [ AGH| , Lemma C.3], |p\ lad . 
Thm. 7.1.9], we present it here with the proof. 

Lemma 2.2 Suppose that M is an open subset of M™ and let K : 
L 2 (M) — > L 2 (M) be an integral operator with the kernel K(-,-). Then 

< \\K\\$ B := I sup / \K(x,x')\dx' sup / \K(x,x')\dx 

\x£M J M x'eM J M 

PROOF: The claim follows from the inequality 

\m P , P <\\K\\{/?\\K\\t/%, (2.7) 
where K is now an integral operator on L P (M), p^ 1 + q^ 1 = 1, and 

||-K"||oo,oo := SU P / \K(x, x')\ dx' , \\K\\x t i := sup / \K{x,x')\dx. 

x&M J M x'eM J M 
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If K is bounded for p = 1, oo, we can prove ( |2.7| ) for the other p by 
an interpolation argument adapted from the discrete case ||Mad|| . By 
Holder's inequality 



K(x, x')ip(x')dx' 



< / \K(x,x')\p\K(x,x')\«\il;(x')\dx' 
'm 



< 



\K(x,x')\\^{x')\ p dx' 



M 



\K(x, x') \ dx' . 



M 



so we can easily estimate the 2/-norm of Kif), 



www 



< 



< 



< 



dx 

M 

K\\p/q 



I K\\ p / q 

I II 00,0 

|_Kl|P/9 



K(x, x')i/j(x') dx' 



M 



dx / \K{x, x')\\ip(x')\ p dx' 



M 



M 

K\ 



dx'\tp{x')\ p / dx\K(x,x')\ 



Al 



1,1 



P 

P ' 



which yields the result. | 

Recall that || • ||sh is not a norm and that it simplifies for the 
symmetric kernels, H-K'Hsh = sup.j. eM J M \K(x, x')\ dx'. We are now 
ready to prove the following key result. 

Lemma 2.3 w t— > M(a(w)) is a bounded and analytic operator-valued 
function, which can be continued from {w £ C | Re w < 0} to a region 
that includes w = 0. 



PROOF: As in [[BGRS| , Lemma 2.2], let Hi C L 2 (tt ) be the space 
of L 2 (f2 ) functions of the form (pxu where tp £ L 2 (]R 2 ). Let further V\ 
be the projection onto this subspace, and Vi := I — V\ the projection 



onto its orthogonal complement in L (fio)- Then Rq 



la 



has an analytic continuation into the region {a £ C|« 2 £ C\[3k 2 ,oo)} 
since the lowest point in the spectrum of H Vi \ ViL 2 (tt ) is k 2 — k\. 
This region includes the domain [0, K\) actually considered. To accom- 
modate the extra factors D,C^, we introduce the quadratic form 
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To check boundedness of this form, it is therefore sufficient to ver- 
ify that R^(a)iViD* and Rq(q)^ViCI are bounded operators. We 
shall check it for their adjoints. To this purpose, it is enough to show 
that C\Vi and DVi are (Rq (a) - 2 "Pj 1 ) -bounded, i.e., that there exist 
positive a, b such that 

W e Q{H X ) : \\C x Vtn < a\\R^{a)-^VM + 6||^|| , 

and similarly for DV±. However, 

iivp^ii 2 = \m+i)^vin 2 -\\nn 2 

\\(H + 1)^|| < ||(g - c?)*vm + VT+^ 11^|| 
< \\R^{a)-^Vi^\\ + VTT^I 



Here V means the gradient in the variables (x,u)) through which 
all the actions of C\,D can be estimated, e.g., \(C\ip)i\ = \A\ijj\ < 
2\\v\\oo\ VV'I, etc. In the same way, one verifies the analyticity of the 
operator- valued function DR^(a)C\, which is equivalent to the ana- 
lyticity of the complex- valued function a. \— > b a {-, •). 

Consider next the regular part of Rq(ol)Vi containing the opera- 
tor N a . Let h be a C°°-function of compact support in M 2 . As pointed 
out above, using integration by parts and the explicit form of the op- 
erators C\, D one sees that it is sufficient to check the boundedness 
and analyticity of hn a h and hn a ^h, where 

n a (x,x') := —K Q (k\(a)\x — x'\) + lnk^a) , 

n a ^(x,x') = : — ki{a)K x {ki{a)\x - x'\); 

2tt \x — x | 

recall that tfi means the derivative w.r.t. x M and K' = —K\ holds true 
- cf. [ |AS| , 9.6.27]. We will use the following estimates which are valid 
for the Macdonald functions ||AS| , 9.6-7] with any z G (0, oo): 



\(K (z) + ln^)e- 2 | < ci , \K x (z) -z- l \< c 2 
\[K 1 (z)-z(K (z)+K 2 (z))/2}\ <c 3 , \zK x (z)\<\. 
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Passing to the polar coordinates, 

X M — x' M = (pCOS<y9,pshl<y9) , p m '■= SUp \x — x'\ , 

a;,a;'Gsupp h 

we check the finiteness of the Schur- Holmgren bounds: 
II^q^IIsh = sup \h(x)\ / \mi(x, x'\ a)h(x')\ dx' 

rpm rpm 

< ciWhWl, / e klia ~> p pdp+ / \\np\pdp 

Jo Jo 

< c 1 \\h\\ 2 00 p m (p m e Klpm + max{e~ 1 ,p m lnp m }) 



\hn a)fl h\\ SR < \\hWlo I = WHlo Pit 



pm pdp 
o P 

Concerning the analyticity, one should investigate the complex-valued 
functions w t— > \6,hn a ( v} )h , ip) and w t— > (</>, hn a ( w )^h ip\ , where <p,ip 
are arbitrary vectors of L 2 (Vl ). Using the Schwarz inequality, it is 
sufficient to check the finiteness of norms of the complex derivative 
w.r.t. w of the corresponding operator-valued functions. Since K[ = 
-(K + K 2 )/2 by [f£SJ, 9.6.29] and h(a(w)) = e w ~\ we put z : = 
ki(a(w))\x — x'\ and write 



dlT'a(w) 

dw 



M = ~(k 1 (z)--) 

Z7T W 2 \ Z J 



dw 2n \x — x'\ w 2 



2 



Using now the inequality w 2 e w < C4 for w G (—00, 0), we are able 
to estimate the Schur-Holmgren bounds: 



dn. 



h- 



a(w) 



dw 



< C2c4 h \\L Pm , 



SH 



dn, 



h- 



dw 



< C3C 4 ||/i||LPm- 



SH 



Thus the derivatives are bounded for w G (— 00, 0), and since the 
limits as w tends to zero make sense, we can continue the function 
analytically to w = 0. | 
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Now we are in position to follow the standard Birman-Schwinger 
scheme to derive the weak-coupling expansion. Eigenvalues of H\ cor- 
respond to singularities of the operator-valued function (J + K"] 
which we can express as 



"a\ — 1 



J+(/ + M A )-% (J + M 



-i 



(2.8) 



Owing to Lemma |2.3| , ||Mx|| is finite and we can choose A sufficiently 
small to have ||Ma|| < 1; then the second term at the r.h.s. of ( |2.8| ) is 
a bounded operator. On the other hand, (/ + M\)~ 1 L\ is a rank-one 
operator of the form (ip, -)ip, where 

A 



-— lnfci(a) xi(u)C* x , 
In 



ip(x,u) : 
<p(x,u) := (I + M x )- l D X i (x,u) 



so it has just one eigenvalue which is 



[ij;,<p) = -— In An (a) 
2tc 




Xi(u)C* x (I + M x )~ 1 Dxi (x,u)dxdu. 



Putting it equal —1 we get an implicit equation, F(X,w) = 0, with 

F(X, w) := w - A / V Xi (u) CI \{I + Mxy'D X i] (x, u) dx du , 

(2.9) 

where M\ has to be understood as a function both of A and w. Ex- 
panding (/ + Ma) -1 into the Neumann series we find 



F w (0, 0) = 1 + , F A (0, 0) = - — (xi, C* D xi) , 



and by Lemma [2.3| we know that F(\, w) is jointly analytic in A, w. 
In view of the implicit function theorem w = w(X) is then an analytic 
function and we can compute the first term in its Taylor expansion: 



> 



^,a(0,0) 1 

~ ^-(Xi^o^Xi) • 



F tW (0,0) 2tt 
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But (C ) n = for n = 4, . . . , 7, B 3 xi = 0, and (A 2 Xi, B 2 Xi) = since 
[L 2 Aw = 0. It follows that 

-jr(0) = 77- (^1X1, = ~- [ Xi(u) 2 du [ v(x) dx = -— (v), 

dX 2ti 71 J J R 2 71 

(2.10) 

where we have employed the symbol (v) := J r2 v(x) dx. 

We note that a 2 — > k\— holds as A — > 0+, and consequently, 
ki(a) — > 0+. Thus w(0) = is well defined because w = (ln/ci(a)) _1 
by definition. Furthermore, the solution a 2 clearly represents an eigen- 
value if and only if w is strictly negative for A small. A sufficient con- 
dition for that is that the first term of the expansion of w(X) is strictly 
negative; due to (|2.10|) it happens if (v) is strictly positive. Summing 



up the discussion, we get the announced three-dimensional analogue 
to Theorem 1.2 in 1BURSJ: 



Theorem 2.4 Let fi A be given by fl2l[), where v G C£°(R 2 ) sat- 
isfies (v ) > 0. Then for all sufficiently small positive X, — A^ A has a 
unique eigenvalue E(X) in [0, nf), which is simple and can be expressed 
as E(X) = k\ — e 2w( - x ^ , where X 1— > w(X) is an analytic function. 
Moreover, the following asymptotic expansion is valid: 



W (X) = -X-^(v) + 0(X 2 ) 



7T 



3 An alternative method 

Now we will derive the weak-coupling expansion by constructing the 
asymptotics for singularities in a particular boundary value problem. 
This approach enables us to derive easily higher terms of the expan- 
sion. At the same time it allows a unified treatment for different di- 
mensions; in this way we will be able to amend the existing results 
concerning deformed strips. 

First we introduce a unifying notation. Let n = 2, 3 be the di- 
mension of the considered deformed region, i.e., the perturbed pla- 
nar strip or layer, respectively. We set x = (x 1 , . . . , x n ~ l ) G M n_1 
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and (x,u) G Q$ := M" -1 x(0, d) for the unperturbed domain. From 
technical reasons it is convenient to change the setting slightly, in 
comparison with (|2.1 ) and ||BGRS|| , ||EV3|| , and to deform the "lower" 



boundary of Qq what we certainly can do without loss of generality. 
We denote therefore in this section 

Q x := {(x, !i)£R": -Xdv(x) < u < d} (3.1) 

with v G C|f ) (]R"~ 1 ). We denote by —A' the (n — l)-dimensional Lapla- 
cian, while —A stands for the n-dimensional one. We also use 



(/>:= / f(x)dx. 



|| ■ || as the norm in L 2 (M n and 

( \ f 171 at+\ ft if n = 2 

aim) := < , , x fi(t) := < . .„ 

[ (mm) [ Int if n = 3 

3.1 The asymptotic expansion 

Let us now construct the asymptotics of the eigenvalues m\ of the 
following boundary value problem: 

(A + n\)^ x = ml^x in fi A 
(x, Xdv(x)) = ^\(x, d) = 

as they approach zero. We will seek it in the form 

IXl if n = 2 

{ ( -A (3 ' 2) 

ex P ( - Xrrii) J if n = 3 



where the existence of such expansions follows from ||BGRS|| and The- 



orem |2.4| , respectively. Notice that this corresponds to the expansion 
of E(\) — k\ — m 2 x , the ground-state eigenvalue of — A^> A in the prob- 
lem discussed above, because the mirror transformation of Q\ on ( |3.1D 
does not affect the spectral properties. 
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Suppose that a function / G C^°(lR n_1 ), supp / fl supp t> = 0, and 
(/) 7^ is given. If we manage to construct a solution if)\(x,u;m) of 
the boundary value problem 

(A + kD^x = m 2 i) X + (a(m) - a(m x )) f\i in fi A (3.3) 

ipx = on d£l\ 

which is bounded and non-vanishing w.r.t. m for small nonzero m, 
then *f!\(x, u) = ip\(x, u; nix). We shall look for the asymptotics of ipx 
in the following form, 



ipx(x, u;m) = y~] X l ipi(x, u\ 



m) 



(3.4) 



i=0 



Substituting ( |3.4j ) and ( |3.2| ) into (|3.3|), we obtain a family of the bound- 
ary value problems: 

(A + kI)i/j = m 2 ip + a(m)fxi in fi i = (3.5) 

■^0 = on dQ 



(A + ^)^ = m 2 ^ + (-l)"- 1 m i /xi in *>1 (3.6) 



ipi — if u = d 
<l J i ryiPc, j 



3=1 



if u = 



One can check easily that ipo — ~ ot{m)(— A' + m 2 ) solves ( |3.5| 

and has the asymptotics 



ip (x,u;m) 



2tt 



(f) 



+ (-l) n - 1 a(m) 



+0 (a(m) 2 ) 



P(\x-x'\)f(x') dx' + 5 6 n (j-ln2)(f) 



(3.7) 



as m — > 0, where 7 is the Euler number and 5 J n the Kronecker delta. 
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Lemma 3.1 Suppose that F e C°°(f2o) with a bounded support and 
H e C^(]R n_1 ) have the expansions 

oo oo 

F(x, u; m) = a{m) l Fi(x, u) , H(x; m) = a{m) 1 Hi{x) 

i=0 i=0 

as m — > 0. Define F^ := , u) Xk(u) du . Let O be the solution 

of the boundary value problem 

(A + /^)0 o = F o in Q , (3.8) 
0o = if u = d, 
0o = H Q if u = 0; 

then the condition 

<F 0) i) = y!«i<ifo) (3.9) 

necessary and sufficient for existence of a solution of the boundary 
value problem 

(A + K\)(j> = m 2 (f) + F in Q , 

= if u = d , 
= # if u = 0, 

which is bounded as m — > 0. If it is satisfied, the solution has the 
asymptotics 



0(x, u; m) 



\n-l 



Proof: The statement is obvious if H = 0. In particular, the solu- 
tion is constructed by the Fourier method in the explicit form 



0(x, it; m) = ^2 <M x ; m )Xi(u). 



i=i 
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By a direct calculation it is easy to see that 0j are bounded functions 
for ra > so long as i > 2. The problem arises for i = 1, because in 
general 0i tends to infinity as ra — > 0. The condition ( |3.9| ) guarantees 
that the explicit solution has no such pole. This proves the suffi- 
ciency. To see that the condition is necessary at the same time, one 
integrates by parts in the scalar product equation 

(xi,(A + «?-m 2 )0) = {xi,F) 

and puts m = afterwards. In the opposite case, H ^ 0, we use the 
replacement 

<f>(x, u; ra) = (p(x, u; ra) + (\ — — j H(x; ra) 

and expand the r.h.s. of the equation for ip in the Fourier series, which 
reduces the task to the previous situation. | 

Corollary 3.2 e C°°(Q) holds for any bounded domain Q C fi . 



It follows from Lemma that the recursive system of the boundary 
value problem ( |3.6|) has solutions which are continuous with respect to 
ra in the vicinity of ra = and decay as |x| — ► oo for ra > 0, provided 
the mi's satisfy the following recursive relations: 



ra; 




3=1 



dv? 



0; 0) 



(3.10) 



In particular, owing to 



and Lemma |3JJ we get 



K 



rai 



n-2 



7T 



(3.11) 



which agrees with the leading term obtained by the Birman-Schwinger 
method in the previous section - cf. Theorem \2.4\ and ( |3.2| ) - as well 
as with the corresponding result fll.lj) in the strip case. 
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3.2 The next-to-leading order 



Let us now calculate m 2 . By virtue of ( |3.6| ), ( |3.7| ) and ( p. 11] ) the 
boundary value problem for ipi together with the boundary condition 
for ip2(x, u; 0) look as follows 



7T 



n-2 



(3.12) 



4>t 



if u = d 



i A d ^ -f n 

Wi = at> — — if u = 

aw 

■^2 = if u = d 



1p2 = 



9m 



if m = d, m = 



with 



9m 



m-l 



x, 0; m) 



2tt 



n-2 



(3.13) 



where 13(f) is the square bracket from (3/7). Hence 

_ a 12 Kid I dipi 



m 2 



-1 



d (f) \ du 



0;0) 



(3.14) 



and it is sufficient to find With eq. (|3.12|) and Lemma [3~T| in mind, 
we consider the following boundary value problem 



(A + = (-1} 



7T 



n-2 



6n = if u = d 



(3.15) 



and seek 0o in the form 
(-1)"- 1 ; 2 



<j>o(x,u) 



2tt«- 2 V d Kl 



■it 



1 ~ d ' ^ S > dv ^~^ x ' u ^ 



(3.16) 
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substituting it into ( |3.15| ), we arrive at the boundary value problem 

(A + K 2 1 ) ip = -d(f){l-^)(A' + K 2 1 )v + 2K 1 ^(v}fxi in fi 

(p = on dQ . 



The Fourier method gives 



fc=2 



^d^K^v + KK-ATH^f-if) 

a Ki 



Lemma |3TT1 an relations (|3.12|) , (|3.13| ), ( |3.15|) , and (|3.16| ) together with 
the last result imply that 



du 



^,0;0) 



•l) n [2 



2tt 



n-2 



d 



X 



7T 



n-2 



TjTT — 1 v TO ri — 1 



/3(|a; — a/|) /(V) dxdx 1 

-1 

+ </> / P(\x-x'\)v(x')dx' - (v) [ P{\x-x'\)f(x')dx' 

+(f) 



3 „(*) + 2 £ [(-A' + k\ - «?) _1 (-A' - k\) v] ( 5 



fc=2 



+5^ ( 7 -In 2) <„) 

7T 



where we have employed also the implication 

-1 



(F) — (— A') F 



2tt 



n-2 
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Substituting this into (|3.14|) we get the sought coefficient: 



m 2 = 3 ( w 2 ) + _I_ J v (x)(3(\x-x'\)v(x')dxdx' 

R"- 1 XK™- 1 
/ 00 \ 



k=2 

+<^( 7 - In 2) (3.17) 

7T 



3.3 The critical case 

As we have pointed out in the introduction, the above result is most 
interesting in the critical case, (v) = 0, when the first coefficient fl3.11| ) 



equals zero and m2 given by ([3.171) determines the leading order. In 
this situation we have the following result. 

Theorem 3.3 Let V G C^°(lR n ~ 1 ) be an arbitrary function such that 
(V) = and 

u(ar)=v(-), ct>0. 
Then the following inequalities hold, 



n n-2 



^||y||2 + ^_||y|||| A V||-2^ (T 2|| V /( A /)-iy||2 > \ 

y 2 Q.K^CT J 



<m2< _>±^_ ( 3|| y || 2 _ 2/c?o- 2 ||V / (A / )- 1 V" 2 



PROOF: In the first place, note that (V) = implies 

/ V(x)(3(\x-x'\)V(x')dxdx' = \\V(A')- 1 V\\ 2 >0 

T n-2 J 



in—i v hp"— 1 



because A'/3(\x\) = 27i n 2 S(x) holds in the sense of distribution. Under 
the stated assumptions, the formula (3- 171) yields therefore 

2 n-1 f 

m 2 = I^WVf - 2kW\\V{A)- 1 V\\ 2 + 2A(a) 
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where 



fc=2 



and it suffices to find suitable bounds on A(o). 

Since the Fourier transformation together with the Plancherel the- 
orem give the estimate 

IIFII 



(-A> + (4-4)a*y l F 



< 



'A - 4)° 2 



(3.18) 



we obtain the upper bound 



A(a)<- \\V\\ 2 + 



2 2 



IVIIIIAVI 



where the numerical factor comes from YlT=2(k 2 ~~ -0 1 
On the other hand, denoting 



U k (x;a) := [(-A' + (4 — 4) ct 2 ) -1 V 



we see that 



V (-A + (4 - k\) a 2 )' 1 (-A - k\o 2 ) V 

(-A + (4 - 4) a 2 ) U k (-A' - 4a 2 )U k 



Integrating the r.h.s. by parts and using (j3.18 ), we get the lower 
bound 

oo 

A(a) = J2 {\\^'Uk\\ 2 + 4(k 2 ~ 2)a 2 ||VU|| 2 - 4(k 2 - l)^\\U k \\ 2 ) 

k=2 

oo oo _ 

>-Y,4(k 2 -iy\\u t \\ 2 >-\wrY, 1 ^r l = - 



k=2 k=2 

which concludes the proof. 

This theorem confirms the spectral picture we got from ( |1.2|) and (|L 
More specifically, m 2 > as o — ► oo so the critical weakly bound state 
exists for sufficiently smeared deformations, and vice versa. In contrast 
to and ( |1.3| ), however, we are able now to tell from ( |3.17| ) for any 
given zero- mean v the sign of m%. 
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